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On the importance of high redshift intergalactic voids 



oo 

o 
o 

(N 
X) 

O 

<N 

S: 

9^ 
6 



(N 
> 

o 



o 

oo 

o 



Matteo Vieli'2, Jorg M. Colberg^'^ & T.-S. Kim^ 

^ INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, 1-34131 Trieste, Italy 
2 INFN/National Institute for Nuclear Physics, Via Valeria 2, 1-34127 Trieste, Italy 
^ CMU, Physics Department, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA 

UMass, Department of Astronomy, 710 North Pleasant Street, Amherst, MA 01003, USA 
^ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 



20 February 2008 



X 
S3 



ABSTRACT 

We investigate the properties of one-dimensional flux "voids" (connected regions in 
the flux distribution above the mean flux level) by comparing hydrodynamical simu- 
lations of large cosmological volumes with a set of observed high-resolution spectra 
at z ~ 2. After addressing the effects of box size and resolution, we study how the 
void distribution changes when the most significant cosmological and astrophysical 
parameters are varied. We find that the void distribution in the flux is in excellent 
agreement with predictions of the standard ACDM cosmology, which also fits other 
flux statistics remarkably well. We then model the relation between flux voids and 
the corresponding one-dimensional gas density field along the line-of-sight and make 
a preliminary attempt to connect the one-dimensional properties of the gas density 
field to the three-dimensional dark matter distribution at the same redshift. This pro- 
vides a framework that allows statistical interpretations of the void population at high 
redshift using observed quasar spectra, and eventually it will enable linking the void 
properties of the high-redshift universe with those at lower redshifts, which are better 
known. 

Key words: Cosmology: observations - cosmology: theory - cosmic microwave back- 
ground, cosmological parameters - quasars: absorption lines 



1 INTRODUCTION 



Over the past decade, properties of voids have been 
more widely investigated, using different observational 
probes an d tracers (mainly i n the l ocal universe , see fo r 
example iHovle fc Vogeievf (|2004|), iRoias et all (| 2004l) 
iGoldberg et all (l2005t).~lRqias et all (I2005D. iHovle et all 



Il2005h. ICeccarelli et al.1 ||2006|). iPatiri et all (l2006al ). 
iPatiri et al.l (|2006bl ), iTikhonov fc KarachentsevI l|2006l )). 
and theoretical - analytical or numerical - models in the 
framework of the A Cold Dark M at ter (ACDM) concordance 
cosmo lo gv (e.g. [Peebles! I^2001^ lArbabi-Bidgoli fc Miiileil 
ll2002h. iMathis fc Whitl (|2002|) iBenson et al.l (|2003|) 
ICottlober et al.l (120031) . Isheth fc van de Wevgaert 



2004 ) 
(200£) 



Goldberg fc Vogelevl ||2004) 



Colberg et al] (|2005|) 



Padilla et al 



Boleiko et al 
2005h 



iFurlanetto fc PiranI ||2 00?) . 'Hoeft et al.' ('200di) . Ibee fc ParkI 



ll2006l), IPatiri et al.1 (2006c). Shandarin et af 

iD'Aloisio fc Furlanettol ||2007| ). ITuIIvI ll2007h. iPark fc Lee 



2007 ). iBrunino et al.l 



20061) 



2003) , iNevrinckl (|2007l ) . iPeebled (|2007l )) 



van de Wevgaert fc Schaap 



The emerging picture is encouraging. On the one hand. 



some results appear to be somewhat hard to understand, 
such as the fact that galaxies observed at the edges of 
and in voids appear to be a fair sample of th e whole 
galaxy population (the so-called 'void phenomenon' . IPeeblesI 
(|2001 ))■ the fact that dwarf galaxies are not found in voids 
contrary to ex pectations from numerical simulations (e.g. 
iPeeblesI l|2007h . or the cold spot obs erved in the Cosmi c 
Microwave Backgroun d (CMB) data (jRudnick et al.1 (|20o3), 
iNaselsky et al.l l|2007h ). On the other hand, other observa- 
tional results, such as those from the 2dF or SDSS galaxy 
redshift surveys, from s tudies of voids in qu a sar sp ectra, 
or the CMB spectrum (|Caldwell fc Stebbinj l|2007l )). are 
in reasonably good agreement with theoretical predictions 
from numerical simulations of the standard ACDM cosmol- 
ogy. Linking the low-redshift properties of voids to those at 
higher redshifts offers the opportunity to constrain the void 
population over a significant fraction of cosmic time and to 
possibly find out more about voids, to ultimately understand 
their role in the context of cosmic structure formation. 

Unfortunately, there are very few observables that con- 
strain the population of voids at high redshift (see, how- 
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ever, iD'Aloisio fc Furlanettj |20o3)). Here, we will use the 
Lyman-Q forest as a tracer of voi ds in the hlRh- r edshi ft 
universe. The Lyman-a forest (e.g. iBi fc Davidsen' (l997|)) 
has been shown to arise from the neutral hydrogen embed- 
ded in the mildly non-linear density fluctuations around 
mean density in the Intergalactic Medium (IGM), which 
faithfully traces the underlying dark-matter distribution at 
scales above the Jeans length. It is a powerful cosmological 
tool in the sense that it probes the dark matter power spec- 
trum in a range of scales (1 to 80/i~^Mpc comoving) and 
redshifts {z =2-5.5) not probed by other observables. How- 
ever, the observable is not the matter distribution itself but 
the flux: a one-dimensional quantity sensitive not only to 
cosmological but also to the astrophysical parameters that 
describe the IGM. 

The so-called fluctuatin g Gunn-Peterson approxima- 
tion iGunn fc PetersonI [l965h relates the observed flux to 
the IGM overdensity in a simple manner once the effects of 
peculiar velocities, thermal broadening, and noise properties 
are neglected: 



F(x, z) = exp[-A(2)(l -I- 57gm(x, z)f] . 



(1) 



with l3{z) ~ 2 — 0.7(7(2) — 1), 7 the parameter describ- 
ing the power-law ('equation of state') relation of the IGM, 
T = 70(2) ((1 -I- SigmY''^ , and A{z) a constant of order 
unity which depends on the mean flux level, on the Ultra 
Violet Background (UVB) and on the IGM temperature at 
mean density. In practice, linking the flux to the density of 
the dark (or gaseous) matter requires the use of accurate 
hydrodynamical simulations that incorporate the relevant 
physical and dynamical processes. However, for approximate 
calculations equation (1), in conjunction w ith simple phy s- 
ical prescriptions of Lyman-a clouds (e.g. ISchavd l|200ll) ). 
can offer precious insights. 

Provided the flux-density relation is accurately mod- 
elled, using the Lyman-a forest to check the properties of 
the population of voids in the matter distribution would 
avoid having to know the bias between galaxies and dark 
matter and would allow to sample over larger volumes the 
values of cosmological parameters that affect the void pop- 
ulation. On the other hand, the fact that the flux infor- 
mation is one-dimensional requires non-trivial efforts to 
relate it to the three-dimensional (dark) matter density 
field. Motivated by huge amounts of new data of exquisite 
quality, we here present a first attempt in this direction. 
Note that in the early 1990s, before the new paradigm of 
the Lyman-a forest absorption was proposed, voids from 
the transmitted Lyman-a flux of QSO spectra were anal- 
ysed, usingdiscrete st atistics on the l ines, b y several groups 
(ej 



Crottj 1IT987I'): iDuncan et all l| 19891): lOstriker et al.1 



1 1988t ): iDobrzvcki fc Bechtoldl (Il99ll ): iRauch et all l|l992M . 

These attempts, however, did not reach conclusive results 
due to small-number statistics. We defer to a future pa- 
per for a link between the observationa l properties out- 
lined h ere and those at l o w red shift (e.g. iMathis fc Whitg 
l|2002h : iHovle fc VogelevI l|2004 )') and for a careful inves- 
tigation of the abso r ber-g alaxy relation (e.g. IStocke et al.l 
l|l995h : lMcLin et al.l tOOm . 

This paper is organized as follows. In Section [21 we 
briefly describe the data set used. In Section [3] we present 
the hydrodynamical simulations. Section|4]contains the bulk 



of our analysis and the results, which are summarized in the 
conclusions in Section [5] 



2 THE DATA SET 

We use the set of 18 high signal-to-noise (S/N> 30-50), 
high-resolution {R ~ 45000) quasar spectra taken with 
the VLT/UVES presented in (|Kim et al.ll2007l ). We regard 
these spectra as a state-of-the-art, high-resolution sample 
in the sense that many statistical properties of the trans- 
mitted flux (mean flux decrement, flux probability distribu- 
tion function, and flux power and bispectrum) have been 
accurately investigated and interpreted by several authors 
fro m the same or fr o m sim i lar data sets (se e, for exam- 
ple, 'McDonald ot al' (2000"), 'Kim et al.l (|2004|) . Ijena et al.1 
(2005),. Viel ct al. (2004b), Boh on et al.1 l|2007h '). More im- 
portantly, systematics effects such as metal contamination, 
noise and resolution properties, and continuum fitting er- 
rors have been extensively addressed. From the data set, we 
discard the four highest redshift QSOs, in order to have a 
more homogenous sample at a median < z >— 2.2 and with 
a total redshift path Az = 6.1. The pixels contributing to 
the signal are in the wavelength range 3282-4643 A while 
the wavele ngth ranges for e ach single QSO spectrum are 
reported in iKim et al. I (Hoot'). Both Damped Lyman-a sys- 
tems and metal lines have been removed from this sample. 
The evolution of the effective optical depth derived from 
this data set, TcS — —\n < F >, which - as we will see 
- is the main input needed to constrain and identify voids, 
has been widel y discussed and c ompared with various other 
measurements (|Kim et al.ll2007l ). 



3 THE SIMULATIONS 

We use a set of simulations run with GADGET-2, a paral- 
lel tree Smoothed Particle Hydrodynamics (SPH) code that 
is based on the conservative 'entropy-formulation' of SPH 
(S pring cl 2005.). The simulations cover a cosmological vol- 
ume (with periodic boundary conditions) filled with an equal 
number of dark matter and gas particles. Radiative cooling 
and heating processes were followed for a primordial mix of 
hydrogen and helium. We assume a mean UVB produced b y 
quasars and galaxies as given bv iHaardt fc Madaul l|l996l ). 
with helium heating rates multiplied by a factor 3.3 in or- 
der to better fit observational constraints on the tempera- 
ture evolution of the IGM. This backgr ound naturally give s 
ri2 ~ 1 at the redshifts of interest here (jBolton et al. 120051 ). 



The star formation criterion very simply converts all gas 
particles whose temperature falls below 10^ K and whose 
density contrast is larger than 1000 into (collisionless) star 
particles (it has been shown that the star formation crite- 
rion has a n egligi ble i mpact on fiux s t atisti cs, see for example 
IViel et al.1 (|2004l ') and lBohon et al.1 (|2007 | ')'). More detaUs of 
the simulations can be found in Iviel et al.ll2004h . 

The cosmological model corresponds to the 'fiducial' 
ACDM Universe with fin, = 0.26, = 0.74, fib = 0.0463, 
Us = 0.95, and Hp = 72 km s~^ Mpc"^ and as = 0.85 (the 
B2 series of IViel et all (|2004f )'). We use 2 x 400^ dark matter 
and gas particles in a volume of size 120 h-^ Mpc box. For 
cross checks we also analyse smaller boxes of size 60, 30, 
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Figure 1. Slice of the projected gas density (thickness is 10 WLpc/h) from our reference ACDM hydrodynamical simulation at z = 2.2. 
The presence of roughly spherical regions below the mean density of sizes of the order tens of Mpc is clear, denser environments at 
S ~ 100 tend to line up in filaments or could be found (more rarely) inside underdense regions (voids). A large underdense region is 
shown as a white thick circles, while haloes with total masses above 2 X 10^^ Mq, that trace the cosmic web, are represented by smaller 
white circles with radii proportional to their masses. The dashed line-of-sight will be analysed more closely in Section [4.4l Filled triangles 
indicate some underdense regions of different sizes: usually smaller voids are surrounded by many massive haloes, while larger ones have 
less haloes at their edges (with smaller masses). 



and 15 h ^ Mpc. The gravitational softening was set to 
in comoving units for all particles for the largest 
simulated volume and changes accordingly for the smaller 
volumes. In the following, the different simulations will be 
referred to by the tuple (comoving box size in Mpc/^,(no. 
gas particles)^''^), so (120,400) denotes our fiducial simula- 
tion etc. (all the sizes reported in the rest of the paper will be 
in comoving units). In Figure [T] we show a two-dimensional 
slice of thickness 10 comoving Mpc/Zi of the gas density for 
our reference simulation at z = 2.2. Large regions whose 
gas density is below the mean are present and overdense 



regions separate them and give rise to the cosmic web. Oc- 
casionally denser regions of overdensity ~ 100 reside in the 
inner regions of voids, however their sizes are quite small and 
their impact parameter in the lines-of-sight to distant QSOs 
should be also relatively small. In Figure [T] we also overplot 
the haloes whose total mass is above 2 x W^^ Mq (white 
circles) and show some underdense regions of different sizes 
(filled triangles) . Usually, the larger voids are surrounded by 
few and not very massive haloes that line up in filaments, 
while smaller voids, such as the one at (x — 77, y = 73) 
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Figure 2. Box size and resolution effects at ^ = 2 on the frac- 
tion of one— dimensional voids in tfie flux distribution (of size 
R or larger), extracted from four different simulation volumes 
(120, 60, 30 and 15 /i~^Mpc) and three different resolutions 
(200^, 400^, 512^ gas particles). The flux has been smoothed with 
a ID top— hat filter of size 1 Mpc, which roughly corresponds 
to the Jeans length at this redshift. Spectra have been normalized 
to the same effective optical depth. 



Mpc//i have more spherical shapes and have more massive 
haloes around them. 

Note that some numerical efforts in order to simulate 
the intergalactic voids are required: on one side large vol- 
umes need to be simulated in order to sample the void dis- 
tribution at large sizes, on the other side relatively high 
resolution is mandatory. In fact, with poor resolution the 
denser regions at the edge of voids are not properly resolved 
and will result in a smaller amount of neutral hydrogen and 
thereby less absorption: this will cause an overestimate of 
the void sizes. 

The parameters chosen here, including the thermal his- 
tory of the IGM, are in reasonably good agreement with 
observational constraints and with recent results from CMB 
data, weak lensing a nd other data set s obtained from the 
Lyma n- g forest (e.g. IViel et al.l (|2006l ). iLesgourgues et al] 
(l2007|)). The mock QSO spectra extracted from the simula- 
tions have the same resolution and noise properties as the 
observed ones, even though this will not affect our results. 
Here, for two reasons we consider only simulation snapshots 
at 2; ~ 2. First, the high-resolution data set has a me- 
dian redshift < z >— 2.2. Second, at these redshifts the 
UVB should be uniform, and the IGM's equation-of-state 
is reasonably measured, leaving little room for ionization 
voids to be labelled as matter voids or for temperature fluc- 
tuations to strong ly affect the Lyman-Q forest properties 
jShang et al.ll2007). However, we caution that: i) recent re- 
sults from Bolton et al.l (|2007l ). based on the flux probabil- 



Figure 3. Fraction of ID voids in the flux distribution at z = 2.2 
of size larger than R as extracted from the observed sample 
(Icr Poissonian error bars). Continuous lines represent the refer- 
ence hydrodynamical simulation (120Mpc//i) with voids identi- 
fied taking ±3cr values around the observed effective optical depth 
value (dashed line). The flux has been smoothed with a ID top- 
hat filter of size IMpc/h. 



ity distribution function, seem to suggest that the thermal 
state of the IGM at 2 < 3 could be more complex than 
expected; ii) patchy Hell reionization at z '-^ 3.2 c o uld a f- 
fect the uniformity of the U VB (e.g. iBolton et al] (|2006h : 
iFaucher-Giguere et al.l (|2007l )) and some effects at smaller 
redshifts might influence Lyman-a opacity. 



4 RESULTS 

4.1 The ID flux void distribution 

We search for flux voids in the simulated spectra in the fol- 
lowing simple way: i) we extract 1000 lines-of-sight from the 
simulated volume; ii) we normalize the mock spectra to have 
the observed mean flux value; Hi) we smooth the flux with a 
top-hat filter of size 1 Mpc, which roughly corresponds 
to the comoving Jeans length; iv) we identify connected flux 
regions above the mean flux level at the redshift of interest. 
The normalization is performed by multiplying the whole 
array of 1000 simulated optical depths arrays by a given 
constant (found iteratively after having extracted the spec- 
tra) in order to reproduce the wanted Toff value. The choice 
of a different mean flux will of course impact on the void 
selection criterion, since a smaller (larger) mean flux level 
will increase (decrease) the number of pixels at high (low) 
transmissivity that are likely to correspond, in a non-trivial 
way, to under(over)-dense regions. 

As a first preliminary test we check for the effects of res- 
olution on the flux fractions of voids at z = 2. Results are 
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shown in Fig. [2] which shows simulations with different box 
sizes and resolutions, but with the same initial conditions 
normalized to the same mean flux. Note that this analysis 
is performed in redshift space and peculiar velocities affect 
the flux properties. From the Figure, it is clear that in order 
to sample the fraction of voids accurately a large enough 
volume and relatively high resolution are needed. For exam- 
ple, the (120, 200) simulation differs dramatically from the 
(120, 400) one, because in the latter relatively dense regions 
placed at the edge of void regions are better resolved. How- 
ever, this discrepancy disappears once the flux is smoothed 
over a given scale (the Jeans length in our case). In a way, 
smoothing the flux degrades the resolution of the simulation 
making it less sensitive to the relevant physical processes 
causing absorption at and below the Jeans scale. Also note 
how increasing resolution beyond that in the (120, 400) sim- 
ulation - the (120, 512) case - produce results that are well 
within the statistical error bars that are observed and will be 
shown in Figure |3l Because of this, (120,400) is now taken 
to represent our flducial run that will be used to explore 
other physically meaningful parameters. 

We now perform the same search for the observed QSO 
spectra. The results are shown in Figure |3] which shows the 
distribution of the flux void population extracted from the 
observations (points with Poissonian error bars) and from 
our largest volume simulation. From the Figure it is clear 
that the flducial (120, 400) run is in excellent agreement with 
observations and that flux voids of sizes > 35 Mpc are 
about 1000 times less common than voids of sizes larger than 
few Mpc. We also note a bump at around 20 Mpc/h in 
the voids fraction, but one would need a larger sample of 
QSOs in order to better test its statistical significance. The 
continuous lines show how the void fraction changes when 
the effective optical depth is varied at a ±3o" level around 
the observed value. In particular, the lower and higher ToB 
values (tch = 0.136 and TeS = 0.19, respectively) have been 
chosen in such a way to conservatively e mbrace at a confl - 
dence level of Scr the values obtained in IViel et al. I l|2004l ). 
where the measured value for the effective optical depth was 
Te ff = 0.163 lb 0.009 . The recently determined power-law flt 
of iKim et al.) l|2007l ) r^a = (0.0023 ± 0.0007)(1 + 2)3 65±o.2i 
by using QSO spectra in the range 1.7 < z < 4, is influ- 
enced by the scatter coming from QSO spectra over a large 
range of redshift, and would in principle allow for a larger 
range of effective optical depths. However, by selecting only 
a subsample at a given redshift, smaller statistical errors on 
the effective optical depth can be obtained. 

We stress that the definition of void used here is based 
on the transmitted fiux and thus is very different from stud- 
ies based on the distribution of voids as inferred from the 
galaxy distribution. 

4.2 Exploring the parameter space 

We now discuss the main cosmological and astrophysical pa- 
rameters that can affect the population of voids in the fiux 
distribution. In Figure |3J we plot the ratios of the flux void 
fractions in the flducial and in other simulations, where we 
changed TcS at z — 2.2 within the observational bounds and 
the value of as- As could be expected, the fraction of voids 
changes dramatically with a change in Tcs'- at z ^ 2 the flux 
probability distribution function (pdf) is a steep function 



3.51 

: (120,400)i^ffmin 

3 

\ ' : (120,400)x,ffmax_.. 

_#2.5; (120,400) Og^l 

"S 
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0.0 \ I,, ',,1,1 
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Figure 4. Ratio of the flux void fraction of the (120, 400) sim- 
ulation (erg = 0.85) with runs that use a different value of the 
power— spectrum amplitude (erg = 0.7 and erg = 1, represented by 
the dashed and continuous lines, respectively). Also shown are the 
ratios when the effective optical depth (which is used as a thresh- 
old to define voids in the flux distribution) of the reference run 
is conservatively varied between the observational bounds (±3cr 
values represented by the dot— dashed and dotted lines). Spectra 
have been normalized to the same effective optical depth. 



around the mean flux level. Thus adopting a different void 
selection criterion has a large impact on the void distribu- 
tion. In the following, we normalize the mock QSO sample 
to reproduce the same effective optical depth. Changing erg 
appears to have a smaller effect. A simulation with erg = 1 
results in a faster evolution of cosmic structures and in a 50% 
increase of the void population of size R = 20-35/1"^ Mpc, 
compared with the reference as = 0.85 case. In the sim- 
ulation with a lower value of as, the void sizes are more 
evenly distributed. Note, however, that the trend at the 
largest sizes (45Mpc//i) is somewhat different: here the low 
erg simulation predicts more voids than the higher erg, since 
the smaller amount of power in the former can occasionally 
produce large regions devoid of matter and absorption. We 
also checked the impact of varying other parameters such as 
f^m. Ho and a change in the shape of the linear dark mat- 
ter power spectrum in the initial conditions to account for a 
presence of warm dark matter (WDM) and adding some ex- 
tra power (EP) at intermediate (from Mpc to tens of Mpc) 
scales, that results in an effect opposite to that of WDM. 
The mass chosen for the WDM particle is 0.15 keV. A mass 
that is already ruled out from the fl ux power spectrum of 
the Lyman-Q forest (see IViel et al] l|2007l ) and that pro- 
duces a suppression of power at scales around 0.4 Mpc/h, 
which are marginally sampled in the initial conditions of our 
(120,400) simulation. However, for the qualitative purposes 
of this paper this simulation should show the impact of a 
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Figure 5. Ratio of tlie flux void fraction of the (120, 400) simula- 
tion (iT8 = 0.85) with runs that have some extra power (EP, con- 
tinuous line) and less power (WDM, dashed line) at Mpc scales. 
Also shown are the ratios when the Hubble constant is set to a 
smaller h = 0.45 value (dot-dashed line) and a simulation with a 
larger value of = 0.4 (dotted line). The reference B2 (120,400) 
run has h = 0.72 and Qrn = 0.26. Spectra have been normalized 
to the same effective optical depth. 



Figure 6. Ratio of the flux void fraction for the (120, 400) runs, 
comparing the fiducial run (7 = 1.6 and T = To) with runs that 
have a different value of the IGM temperature at mean density 
(2 To and To are represented by dashed and continuous lines, re- 
spectively) and with a run that has a different value of 7 (dotted 
line). Spectra have been normalized to the same effective optical 
depth. 



WDM particle on the distribution of the largest voids in the 
flux. At the scale of 2.5 Mpc//i the power is 100 times higher 
(smaller) in the EP (WDM) simulation than the correspond- 
ing (120,400). The results are plotted in Figure [5] where all 
the simulations have the same value of as — 0.85: it is clear 
that the WDM and EP runs have opposite effects. WDM 
significantly enhance by a factor larger than 3 the fraction 
of largest voids, by erasing substructures below a given scale, 
while the EP simulation presents more collapsed structures 
that suppress the number of the largest voids. The effect of 
a larger value of fim = 0.4 is similar to that of the EP simu- 
lation and determines a 90% reduction in the fraction voids 
whose sizes are larger than 30 Mpc//i. Finally, we checked 
for the effect of a smaller value of Hp = 45 km /s/Mpc a 
value proposed recently by [Alexander et"al] |20o3) to be the 
true value of the Hubble constant outside a local underdense 
region of very similar density to those explored here. Hav- 
ing a smaller value for the Hubble constant produces effects 
that are very similar to those of a smaller value of as = 0.7 
(see Fig. |4]): the evolution of the cosmic web is less pro- 
nounced and can occasionally produce large regions devoid 
of absorption, even if on average the effects are quite mild. 
We also varied the spectral index but found small dif- 
ferences to the reference case. As can be appreciated from 
Figure [3l the statistical errors on the observed void fraction 
are somewhat large and thereby distinguishing between dif- 
ferent cosmological model is still very difficult and none of 
the models presented here could be convincingly ruled out. 



Several other parameters could have an impact on the 
void distribution, in particular the two parameters that are 
poorly constrained and that describe the thermal state of 
the IGM: To and 7. We investigate the effect of a different 
value of To in Figure |S] A colder (hotter) IGM will result in 
a higher (lower) density of neutral hydrogen and thereby in 
more absorptions, which, in turn, will decrease (increase) the 
fraction of large voids. In Figure|Sl we also overplot the effect 
of changing 7 from the reference value (7 — 1.6) to 7 = 1.1 
(represented by the dotted line in the Figure), which is a 
better fit to the data at this redshift (|Bolton et al.l 120071 ). 
A flatter equation of state will not only make underdense 
regions hotter, but it will also make the slightly overdense 
regions colder. The overall effect is a slight decrease in the 
fraction of flux voids, since cold regions embed on average 
more neutral hydrogen than hot ones. This result shows that 
by selecting mean-flux-level regions we are sampling regions 
around mean density (for which the power-law 'equation of 
state' could be a poor fit to the actual distribution). 

Note that the effect of changing 7 was addressed by 
running a new simulation that self-consistently results in a 
low 7 at 2 = 2 (using a modified version of the GADGET-2 
code kindly provided by J. Bolton) and not by rescaling the gas- 
temperature relation a posteriori as is sometimes done. Another 
important input parameter is the overall amplitude of the UVB, 
paramet rized by r_i2, for w hich we have some observational con- 
straints llBolton et al.ll2005h . However, this quantity is degenerate 
with the value of T^g, so changing r_i2 will result in the same 
change as changing r^ff. 
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Figure 7. ID gas density and temperature in void regions along lincs-of-sight extracted from the fiducial simulation at z = 2.2. Voids 
have been identified from the fiux distribution as regions above the mean fiux level. Results are shown with (triangles) and without 
(diamonds) peculiar velocities, and error bars represent the scatter. 



4.3 From flux to ID gas density 

Having identified voids at z ~ 2, using the mean flux as a thresh- 
old, we are now interested in linking the flux properties to the 
underlying physical properties of the gas density field. The main 
problem we have to deal with is that the flux is observed in red- 
shift space, while we want to recover gas properties in real space. 
The role of peculiar velocities in determining the properties of the 
transmitted flux has been investigated by several authors (e.g. 
IWeinbere et al.l l ll998t) ). however for the purpose of the analysis 
performed here it produces a shift of the absorption feature from 
the real space overdensity to the redshift space observed flux. This 
increases the scatter since redshift space pixels are not related in 
a simple manner to the underlying real space density field. 

Occasionally, peculiar velocities could be of the order of a 
few hun dred km/s, so the c orrespondence is non-trivial (see for 
example I Ranch et al.l l(200^). However, we can rely on hydrody- 
namical simulations in order to infer, at least on a statistical basis, 
the physical properties of the gas. In Figure [7] we show the mean 
values and scatter (r.m.s. values) in the real space gas-overdensity 
and temperature fields, corresponding to the flux voids of the pre- 
vious sections. We show the fields both without (diamond sym- 
bols) and with (triangles) peculiar velocities, which alter the flux 
distribution observed in velocity space. From the plot, it is clear 
that selecting regions above mean flux at z ^ 2, as we did in 
the previous Sections, allows us to identify and measure (with a 
large error bar) the density in ID gas voids with sizes larger than 
~ 10 Mpc in the case with peculiar velocities. If peculiar ve- 
locities are neglected the scatter is much smaller (as expected), 
and the density can be measured for somewhat smaller voids (sizes 
larger than ^ 7 Mpc). For even smaller voids only upper lim- 
its of the order a few times mean density can be set, and very 
little can be inferred about their temperature. The gas tempera- 
ture inside the largest flux voids is of the order of lO"* K. These 
regions are thus cold and, as expected from the nature of Lyman- 
a forest absorptions, sample the IGM 'equation of state' around 
mean density. Therefore, even when peculiar velocities are taken 



into account, the physical properties of the largest voids can be 
reliably studied. Furthermore, note that the effect of the peculiar 
velocities for the recovered mean density and temperature inside 
voids is systematically slightly higher than the corresponding val- 
ues for the case when peculiar velocities are neglected. As before, 
this effect is expected since the peculiar velocity shift can move 
some denser clumps of gas into voids regions that increase their 
overall mean density. We stress that the 'gas voids' defined in 
this Section, as identified from the transmitted flux, do not corre- 
spond to underdense regions in the gas distribution but to regions 
around the mean density. 

In Fig.[8]we plot the voids' densities for voids whose sizes are 
larger than 7h~^ Mpc, with and without peculiar velocities and 
with a varying flux threshold. The scatter is large, but for voids 
regions above the mean flux level in 68% of all cases these regions 
have overdensities in the range 5 = [—0.8, 0.8]. If the flux thresh- 
old is lowered, large voids regions start to sample environments at 
higher overdensities. The dashed line represents a linear fit to the 
points (with peculiar velocities): 5 + 1 = p/ < p >= —0.9F/ < 
F > +1.9. This relation should describe the mean ID density of 
large (> 7 comoving Mpc//i) voids identified in the flux distri- 
bution at z ~ 2, as derived from high-resolution Lyman-a forest 
QSOs spectra. We note that a systematic trend is present and 
the recovered gas density is larger in the case with peculiar ve- 
locities: we attribute this to the fact that with peculiar velocities 
the corresponding gas absorptions are shifted and a void region 
identified in the flux is not simply related to an underdense region 
but can incorporate regions at higher density that will increase 
the recovered global density. 

4.4 From ID gas density to the 3D dark matter 
distribution 

In this Section, we link the ID properties of fiux voids to the 
3D dark matter properties. We use a void— finding algorithm for 
the dark matter distribution of the reference simulation. The al- 
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Figure 8. Densities of voids (larger than 7 Mpc) as a func- 
tion of flux threshold level at 2 = 2.2. Results are shown with 
(triangles) and without (diamonds) peculiar velocities, and error 
bars represent the scatter. The dashed line shows a linear fit to 
the case with peculiar velocities. 



gorithm is described in more detail in IColberg et aP (120051 ). Its 
main features can be summarized as follows. The starting point 
for the void finder is the adaptively smoothed distribution of the 
matter distribution in the simulation. Proto-voids are constructed 
as spheres whose average overdensity is below a threshold value 
St, which is a free parameter. These proto-voids are centered on 
local minima in the density field. Proto-voids are then merged 
according to a set of criteria, which allow for the construction of 
voids that can have any shape, as long as two large regions are 
not connected by a thin tunnel (which would make the final void 
look like a dumbell). The voids thus can have arbitrary shapes, 
but they typically look like lumpy potatoes. For each void, we 
compute the radius Rom of ^ sphere whose volume equals that 
of the void. While strictly speaking Rdm is not the actual radius 
of the void, it provides a useful measure of the size of the void. In 
this work here, we adopt a value of 5t = —0.5 (for a justification of 
this choice see below). As a first, preliminary check, we compared 
the void size distribution of th e fiducial simulation with earlier 
results bv lColbere et akl ||2005| ) and found good agreement. Note 
that the size of the grid used for the density field limits the sizes 
of the smallest voids that can be found with this void finder - in 
our case lh~^ Mpc. However, as discussed above, peculiar veloc- 
ities make studies of voids smaller than 7 Mpc very hard, so 
we need not be concerned about the smallest voids at all. 

As a qualitative example, in Figure |9] we show one linc-of- 
sight through the simulated volume. We plot the ID gas density 
field in real space (top panel), the simulated fiux (middle panel) 
in redshift space, and the gas peculiar velocity field in km/s (bot- 
tom panel) in real space. The x-axis indicates the comoving (real 
space) coordinate along the line-of-sight. Flux absorption is usu- 
ally produced by gas overdensities shifted by the peculiar velocity 
associated to the same gas element. For example, the absorption 
in the flux at 70Mpc/h is produced by the gas density peak at 
around 68Mpc//i, shifted by ~ 200 km/s - which corresponds to 
roughly 2h~^ Mpc. 



The shaded regions show the intersections with the line-of- 
sight of 3D voids in the dark matter distribution identified with 
the void finder algorithm described above. In the middle panel 
we overplot the mean fiux level as a dotted line, with the flux 
smoothed on a scale of 1 Mpc, as in Section l4.1l We find that 
the choice of threshold overdensity 5t = —0.5 for the 3D void 
finder results in a good correspondence of the resulting voids with 
voids defined in the flux distribution. From the bottom panel of 
Figure [9] it is clear that the void regions are expanding, and the 
gas velocity flelds for the largest voids are of the order of ± 200 
km/s (see the void centered on around 20Mpc/h). These high 
redshift voids in the IGM will keep growing, get emptier of mat- 
ter and could be the prog enitors of low er redshifts voids of sizes 
similar to the Tully void iItuUvI HoqI)). The peculiar velocities 
along the line-of-sight rise smoothly from negative values to pos- 
itive ones and are sandwiched by regions in which the peculiar 
velocity shows a negative gradient, which could be the signature 
of a moderate shock. It is worth emphasizing that these velocity 
profiles are not always completely symmetric because the line- 
of-sight does not necessarily pierce the 3D voids in their centers. 
Moreover, from the top panel it is clear that the ID density pro- 
file inside the void is quite complex, showing small density peaks, 
which correspond to small haloes that most likely host the ob- 
served void galaxy population. 

In order to provide a more quantitative picture, in Figure [TOl 
we plot the underlying mean values of ID quantities obtained 
once the 3D voids have been identified (using St = —0.5). The top 
panel represents the probability distribution function of the sizes 
of such ID fiux voids, while in the other three panels each point 
shows, from top to bottom, the mean gas overdensity, flux and 
peculiar velocity field for a given 3D DM void. Note that while in 
Section [4.1l we were interested in measuring the properties of voids 
once the mean fiux threshold was set as a criterion to define voids, 
here we use the 3D DM distribution to look for the corresponding 
ID (Lyman-a forest) fiux— related quantities. 

The overall picture suggests that large voids (> 7 — 
lOMpc/h) are usually related to fiux values above the mean, and 
the corresponding mean values of gas density and peculiar veloci- 
ties have less scatter than smaller voids. The corresponding mean 



value for the gas overdensity is (5| 



GAS 



-0.43, so the 3D void 



region is slightly denser in gas than in dark matter, something 
that can be expected since baryons feel pressure an d are thus 
more diffused than the coUisionless dark matter (e.g. IViel et al] 
(2002)). The different value of ^Q^g than the one found in Sec- 
tion 14.31 does not contradict our earlier findings, since now the 
requirement to define a void is less stringent. Now, the actual 
void has to be a quasi-spherical 3D region in the DM distribu- 
tion and not a connected ID region with flux above the mean 
flux. This means that in ID one is more sensitive to clumps of 
matter producing absorptions, while the same clump will have a 
smaller impact on the mean density of the 3D region surrounding 
it. Cases for which the flux voids correspond to actual DM voids 
are avaraged out statistically by small lumps of gas around the 
mean density along the line-of-sight, that cause absorption and 
do not alter the global 3D DM density of the void region 

For this population of large voids the typical scatter in the 
peculiar velocity fleld is 170 km/s. We thus have shown that the 
population of ID large flux voids, i.e. regions above the mean flux 
level, as estimated from a set of mock high resolution spectra at 
z ^ 2, traces reasonably faithfully a population of 3D dark matter 
voids of similar sizes, with a typical St = —0.5. 

We have also checked that the scatter plots of Figure [T0I are 
not very sensitive to the cosmological model (i.e. different values 
of (Tg), which means that the voids physical properties are the 
same for different values of the power spectrum amplitude. 
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Figure 9. Relation between 3D dark matter voids and pliysical quantities along the line— of— sight to distant QSOs. The intersections of 
3D real-space voids in the dark matter distribution along the line— of— sight are represented by shaded rectangles. The top panel shows 
the ID real— space gas overdensity, the middle one the simulated Lyman-a flux in redshift space, the bottom panel the peculiar velocity 
field in km/s. Absorption features in the flux distribution correspond to overdensity peaks shifted by the peculiar velocity of the bottom 
panel (100 km/s roughly corresponds to 1 Mpc comoving). The dotted line in the flux distribution sets the criterion for the void 
definition, while the dashed one represents the flux smoothed over 1 Mpc. This particular line-of-sight is taken along the x-axis and 
at a y-coordinate of 93 Mpc/ft, in the slice of Figure [T] 



5 CONCLUSIONS AND PERSPECTIVES 

We used Lyman-a forest QSO spectra to constrain the void pop- 
ulation at z ~ 2 at a range of scales and redshifts, which cannot 
be probed by other observables. The main conclusions can be 
summarized as follows: 

• the properties (sizes) of voids identified in simulations from 
the flux distribution as connected regions above the mean flux 
level are in good agreement with observed ones; 

• using a set of hydrodynamical simulations and varying all the 
astrophysical and cosmological parameters, along with checking 
box size and resolution effects, we find that the fiux void size 
distribution is a robust statistics that depends primarily on the 
fiux threshold chosen to define the voids; 

• the flux seems to be a reliable tracer of the underlying dark 
matter distribution, altering the statistical properties of the dark 
matter density field (changing amplitude, shape, adding or sup- 
pressing power in a wavenumber dependent manner) changes the 
properties of the void population although at present is difficult to 
give constraints on cosmological parameters using these statistics; 

• fiux voids, i.e. connected regions above the mean fiux levels, 
correspond to gas densities around the mean density, regardless 
of the role of peculiar velocities, which contribute to a scatter in 
this relation, but do not alter significantly the mean values; 

• linking ID voids to the corresponding 3D dark matter voids 
is more difficult, and we used a void-finding algorithm for that: 
we find that 3D DM voids with a mean overdensity of 5f = —0.5 
correspond to the fiux voids that we have defined from the QSO 
spectra. However, the correspondence is good only for voids with 
sizes larger than about 7 — 10 Mpc/h. 

In this paper, we made a preliminary attempt to link the 
population of voids in the transmitted Lyman-« fiux to the un- 
derlying ID gas density and temperature and 3D dark matter 
density. The use of Lyman-a high— resolution spectra is important 



in the sense that it explores a new regime in scales, redshifts and 
densities which is currently not probed by other observables: the 
scales are of order few to tens of Mpc, the redshift range is between 
z = 2—4, while the densities are around the mean density. Further 
studies, that could possibl y rely on wider dat a sets such as the 
low— resolution SDSS data (Shang ct al. (2007) ), on othe r futur e 
observables at higher redshifts (D'Aloisio & Furlanettol | |2007l) ) 
or on tom ographic studies i nvolv i ng QSO pairs or mu ltiple line- 
of-sights jP'Odorico et al.l l|2006h . ISaitta et al.l l|2007t) 1 would be 
important in understanding the dynamical, thermal and chemical 
evolution of the void population and the interplay between galaxy 
and the IGM over a large fraction of the Hubble time. 
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